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ABSTRACT 

Spectral element methods are p-type weighted res dual techniques for partial differential equa- 
tions that combine the generality of finite element nr ethods with the accuracy of spectral methods. 
We present here a new nonconforming discretizatioa which greatly improves the flexibility of the 
spectral element approach as regards automatic me *h generation and non-propagating local mesh 
refinement. The method is based on the introduct on of an auxiliary ‘‘mortar” trace space, and 
constitutes a new approach to discretization-driven domain decomposition characterized by a clean 
decoupling of the loced, structure-preserving residua' evaluations and the transmission of boundary 
and continuity conditions. The flexibility of the mortar method is illustrated by several noncon- 
forming adaptive Navier-Stokes calculations in complex geometry. 


This research was supported by the National Aero, antics and Space Administration under NASA Con- 
tract Nos. NASl-18107 and NASl-18605 while the first author was in residence at the Insitute for Computer 
Applications in Science and Engineering (ICASE), NA lA Langley Research Center, Hampton. VA 23665. 
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Nonconforming Mortar Element Methods: 

Application to Spectral Discretizations 

1 Introduction 

Spectral element methods [22,25,27] are weighted residual techniques for the approximation of 
partial differential equations that combine the rapid convergence rate of spectral methods [6,14] 
with the generality of finite element techniques [8,12,29]. The spectral element discretization, 
coupled to fast order-independent iterative solvers [21,28,32], yields numerical algorithms which 
have proven computationally efficient on both ser al and parallel processors [10,11]. Although the 
spectral element method is, by construction, applicable in complex geometries [16,18,27], the large 
indestructible geometric unit associated with high-order brick elements leads to a certain lack of 
flexibility as regards automatic mesh generation, adaptive mesh refinement, and the treatment of 
moving boundaries. In this paper we present a new method, the “mortar element method” , which 
largely eliminates this rigidity by allowing for nonconforming matching between subdomains. 

The “mortar element method” represents a new domain decomposition approach [7,13] in which 
there is a clean decoupling of local-structure-pres< rving internal residual evaluations and the trans- 
missions of boundary (or continuity) conditions. The method is not based on Lagrange-multiplier 
interface constraints e.g. [9], but rather on the ex, . licit construction of the appropriate nonconform- 
ing space of approximation through the introduction of a new intermediary mortar trace space. The 
explicit-space approach is more appropriate for fest iterative solution than the Lagrange-multiplier 
methods, as it avoids the necessity of solving a coupled, potentially ill-conditioned problem. Al- 
though we develop the mortar methods here for spectral element discretizations, they are also 
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appropriate in the h-type finite element context [4], in which they constitute an extension and 
generalization of classical nonconforming methods [8,9,29,31], 

We present here the “mortar element method” in its simplest form for the solution of two- 
dimensional second-order elliptic and saddle problems. The emphasis is on the numerical formu- 
lation, implementation, and demonstration of the technique, and the illustration of the fiexibility 
of the nonconforming paradigm; theoretical support for the method is given in [4], in which the 
optimality of the discretization is proven. The outline of the paper is as follows. In Section 2 we 
present the basic discretization for the Poisson equation in terms of the function spaces over which 
the standard variational form is to be tested. In Section 3 we present the associated nonconforming 
bases and the resulting set of discrete equations. Conjugate gradient iterative solution of the mor- 
tar discretization is described, illustrating the strong domain decomposition nature of the residual 
evaluation procedure. In Section 4 the extension of the method to the solution of the Stokes and 
Navier-Stokes problem is presented. Lastly, in Section 5 we give several numerical examples. 


2 Spectral Element Nonconforming “Mortar” Spaces 

2.1 Problem Formulation 

W« consider Brst the solution of . Poisson equation on a domain II of Find n(x, „) such 
that 


-V*u = / inn. 


(la) 


u 


0 on 3n, 


(lb) 
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where 30 is the boundary of O, and / is the prescribed force. We suppose that O is rectangularly 

decomposable, that is, that there exist rectangular fubdomains O*, k = such that 

K 

n=[Jn*, vfc,/, jfc n‘nn' = 0. ( 2 ) 

k=l 

The problem (la, lb) is well posed in ^ in thu sense that the following weak formulation of 

the problem admits only one solution: Find u X 1 uch that 

( Vu, Vv) = < /, V > , Vv G X (3) 

Here (.,*) represents the inner product, and < ., > denotes the duality pairing between X and 
its dual space. For the definition of standard spaces norms and inner products we refer the reader 
to [1]. 

For the Galerkin numerical approximation of problem (la, lb), we test the variational form (3) 
with respect to a family of discrete finite dimensioral spaces Xh, where h denotes a discretization 
parameter: Find Ufi G Xh such that 

{VuH, Vvh) = < /, V,. >, Vvfc G Xh^ (4) 

In the case of a conforming approximation, for whicl X^ C X, the convergence and convergence rate 
of towards u is deternuned essentially by stability (ellipticity and continuity) and approximation 
theory (infimum of l|ti — v^Hi.n over all Vh € X^, whore || • ||i,n refers to the norm over fl). In the 
case of nonconforming approximations, for which Xh ^ X, we must also consider the consistency 
error, which measures the deviation of the approximation space Xh from the proper space X [8,29]. 

To date, spectral element approximations [22] have been based on domain decompositions that 
satisfy (2) as well as the additional constraint th it the intersection of two adjacent elements is 
either an entire edge or a vertex; this second constraint is derived from the conforming assumption, 
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and is also present in the finite element method. In the spectral element context this constraint 
can be prohibitively restrictive due to the large geometric units involved. Although relaxing the 
conforming constraint clearly introduces a new source of error, it has the potential advantage of 
greatly increasing the flexibility of the numerical method as regards mesh generation and adaptive 
refinement procedures. This increase in flexibility improves not only the efficiency of the algorithm, 
but also the tractability of cedculations involving moving and sliding meshes [15]. Furthermore, 
the nonconforming approach achieves generality at no cost in loss of local structure, an important 
consideration as regards optimal solvers. 

We present here a spectral method based on nonconforming approximations in which the con- 
sistency errors are commensurate with the approximation errors. To present the nonconforming 
spectral element space Xh we first describe the anatomy of the discretization. The K rectangular 
subdomains of (2) are now identified as spectral elements, and the (i, y) coordinate system is chosen 
so as to be aligned with the edges of the fl*. These edges are denoted F**', / = 1 ,..., 4 , such that 

4 

= IJr ’ . We next introduce the set of “mortars” 7**, where 
1=1 

7*’ = int(n‘ n if) (5a) 

for some k and /, or 

7'’ = int(fT* n an) (5b) 

for some index m, where p is an arbitrary enumeration p = of all (A:,/) and m such that 

mt(n* n n*) or int(n™ n an) is not empty. The intersection of all closures of all 7P defines a set of 
vertices V composed of all elements 

v« = f"ri7" (6) 

where 9 is an arbitrary enumeration q = of all couples (m,n) for which (^ n ^) is not 
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Figure 1: Subdomains and Mortars ol a Nonconforming Decomposition 

pty. (The set I' is equal to the set of all the vertices of the D* by definition of the mortars), 
tly, we define the skeleton S of the mortar syst< m by 

s = (j 7 P = IJ an*. (7) 

p=l Jk=l 

! geometry of the nonconforming decomposition is shown graphically in Figure 1. 

In order to define the nonconforming space we first require an auxiliary mortar space Wh 

Wh^{4>^c°{S), Vp = i,...,jw, ^ 1 ^, € PjvC'r'’), ^1,0=0} ( 8 ) 

5 re Pjv( 7 ^) is the space of all polynomials on ^ ^ of degree < N, The nonconforming space is 


a given by 
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X^= {v G L^(n), V/c = G Pjv(n*) such that 3<f>EWh for which: 

Vg = VA: = such that is a vertex of fl*, (9a) 

and V/ = l,...,4, VA = l,...,iir, G Piv-2(r*‘0, = ^ >* (9b) 


Here Pjv(n*) denotes the space of all polynomials on of degree < N in each spatial direction; 
the spectral element discretization parameter is the couple h = (if, iV). For a conforming approx- 
imation Xf^ is the standard spectral element space; here, and elsewhere in this paper, we assume 
the reader is familiar with the conforming spectral element method [22] . 

Let us summarize the properties of the approximation space Xf^. First, as regards the uniqueness 
of the solution, we note that the uniqueness of the mortar element <f> G is not of major importance 

£LS long as its image G Xf^ is unique; it is not the mortar element, that must be close to 
ti. The uniqueness of the discrete solution follows from the elliptic! ty of the Laplacian form 
(Vu^, Vu/j), Vu/j G Xhj Vv/i G Xf^ with respect to the following “broken iT^(n) norm”, 


r K 1 


( 10 ) 


Although the proof of elliptic! ty is quite involved (see [4]), an elementary proof of uniqueness can 
be readily derived. To wit, we note that if Uf^ and are two solutions of (4), we get 

0 = (VvA, Vvfc) = ^ f «/! = “fc - 

ib=l''^ 


and thus is piecewise constant. Using the fact that the elements of X^ vanish over dfi and are 
continuous at the vertices of "V, it follows that v/i = 0 and thus u^ = 

Although uniqueness of ^ G is not necessary, it is nevertheless true that spurious (or 
parasitic) modes in (f> correspond to unprofitable work, and can potentially cause problems in the 
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p. — < 

► 

n* 

n* 

0* 


Figure 2: Nonconforming Discre- ization Derived from the Refinement 
of a Conforming Approximation 

subsequent solution of the discrete system (see Section 3.2). There is one situation in which the 
uniqueness of <j> follows easily; this is the case wher i for each 7 P there exists an element D* that 
accepts as an entire edge (see Figure 2). This arises, for instance, from a refinement of a 
mesh which is initially conforming. In this paper ,ve shall consider only this “refinement” case; 
development and analysis of the general problem < f Figure 1 is more involved, and is relegated 
to future publications. For the “refinement” case uniqueness of <f> results from the fact that the 
mortar element 4> coincides exactly with the restriction over 7 ^. To show this we note that, 

by construction, the elements <j> and vj = (ui,)|^* co ncide at the endpoints of 7 ". This implies that 
^ ^ is a polynomial of the local coordinate 5, 

where $ is a polynomial of degree < N - 2. Hen , and in what follows, a = x (or S = y) for a 
horizontal (or vertical) mortar, where x (or y) is a nortar-local variable which scales x (or y) such 
that 7 p corresponds to ] - 1 , 1 [ (similarly, s = x (or 5 = y) for a horizontal (or vertical) edge, where 
X (or y) is an element-local variable on ] - 1, l[ whic h scales x (or y) to the appropriate F^-'). From 
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the orthogonality of ^ — to all elements of Pjy_ 2 (' 7 *’) (9b), it follows that $ is necessarily zero, 
and thus <f> is exactly the trace of one piece of v^. The uniqueness of the solution uy, to problem 
( 4 ) thus yields the uniqueness of the corresponding mortar element. 

Let us consider now the consistency error. The scheme (4) based on the definition (9a, 9b) 
of X/, is optimal in that the consistency error is maintained small by the combination of the 
condition (9b) and the vertex condition (9a). In essence, the condition ensures that the jump in 
functions is small in the interior of internal boundaries, whereas the vertex condition ensures exact 
continuity at cross points where the normal derivative has more than one sense. We note that the 
superiority of the L^- (versus pointwise-) matching of and ^ has been demonstrated previously 
[2]. The mortar methods are different from previously proposed nonconforming approximations 
in that the latter are mortarless master-slave spaces, whereas the current approach is democratic; 
this allows for very simple implementation in arbitrary topologies. 

Lastly, the approximation properties of the speice are similar to those of past nonconforming 

approximations. For example, for the case of a square domain decomposed into several elements, 
as a first result one can use the best global polynomial approximation as a bound for approximation 
errors. The combination of stability, consistency and approximation result in an optimed scheme, 
the details, and degree of locality of which, are described in [4]. We note that for the special case 
of infinitely smooth solutions, u/, approaches u exponentially fast as W 00 for fixed K (spectral 
convergence) . 
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3 Representation and Discrete Equations 


3.1 Bases 

Although the spaces Wh and Xh appear quite complicated, they have a simple basis and evalua- 
tion procedure which yields an efficient domain dec< mposition algorithm. In this section we discuss 
the basis, and in the following section we describe esidual evaluation. 

To begin, we write for the space Wh^ 

VpG{l,...,M} (12) 

j=0 

where we assume that all indices increase with increasing x,y. Here the hf are Lagrangian inter- 


polants defined by 


/if gPk(]-i,i[). = % v.-,je{0,...,lV}= 


(13) 


where the ^,(= ) are the JV -|- 1 Gauss-Lobatto Legendre points defined by the zeroes of 

L^(z)(l - z*), and Ln is the Legendre polynomiEl of order N [30] so that 

I (1 - g ^,y(;) VjG{0,...,A}. 

N{N+l)LNiii) z-i 

The definition (12) is not sufficient given the requirement that ^ G must be C°(S); to indicate 
the continuity condition, we resort to diagrammatic methods. The mortar conventions are described 
in T^ble la, with the basis for shown in Figure 3b for the nonconforming mesh of Figure 3a. 
We next construct a representation for v G X,u in terms of the mortar. To begin, we write 


N N 


-U. = E E « e { 1 . .... K} 

1=0 i=o 


(14) 


where the /if are defined in (13). The internal degrees-of-freedom, t,j G {1, ...,N l}^ are 

clearly free, however the boundary degrees-of-freedom are constrained through (9a, 9b). Based on 



10 


□ 


□= 


=□ 


open or dashed symbol/ 
solid symbol 


^(v«) 

j = 0,...,N 

eissign vertex 

source/ destination 
or degree-of- freedom/slave 


(S) 






□ 

□ 

□ 


Q 


Q 


I 


Q 


■^i 


VOO.VOJV, ViVOjfJVAr 

VOi,VjQ,VNj,Vjtf 

j = 1 

vij t,j e - 1}^ 


assign vertex 


mortar-to-edge projection 


sum vertices 


edge-to-mortar sum 


= Eiioi:5Lo<'SA?'(*)'>}'(y) 


b) 


Table 1; Symbols for Diagrammatic Basis Representation 




[ 

f ff 




I i 

X - - X 




f IT 



1 

1 1 1 


a) b) 

Figure 3: a) Nonconforming Mesh and b) Associated Mortar Basis 
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Figure 4: Diagrammatic Representation of the B;tais for Xh on Nonconforming Mesh of Fig. 3a 

the diagrammatic conventions of Table 1 the admissible v are given by Figure 4, where Q derives 
from the projection (9b). In order to construct Q ve require a basis for V’, which we choose as 


where 

To calculate the projection of (9b) we then perform (here exact) piecewise Gauss-Lobatto quadra- 
ture on N -H 1 points on the element edges and mortar segments, giving 

= vie{i,...,w-i} (17) 

j=i ;=o 

where for the destination edge F*-' (u,) and sourc ; mortar 7 P (<^,) 




(18) 
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7'’ 


p*,J 


So 


Figure 5: Definition of Mortar Offset 


and 


LA’’'’ 


P - — ’ 

^ij — 


I r*-' 


I 7" I 


- < 


if«o=0 

2 V^~'^i^)pNSN,j if ao+ I 7" 1=1 r*>' I 


r*-' L 7V-2 


(19) 


Here is the offset of the mortar 7 P from the edge F*-', as shown in Figure 5 ; the endpoint terms 
of (19) derive from the vertex-pinning condition of ( 9 a). Finally we arrive at 


Vi e - 1}, Vj G { 0 , ( 20 ) 

Note that by proper choice of the basis for tp we can explicitly form the matrix Q, that is, we are 
able to directly invert the diagonal inner product B. The alternating-sign term in assures 

that the entries of B are positive. 

Although in practice we shall evaluate from the diagram without forming the global linear 

projection operator, it is nevertheless of theoretical interest to remark that the diagram is equivalent 
to 
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V* , . , VA: 

*J \int«rior 


y ^ijUnUrface^ j 







1^1 

0 


y 

1 ° 

m j 

\ 


V* I . , Vfc 

*J I interior 


or 


Vg, Vp, VjG{l,...,iV-l} 


V * = Q 


(21a) 


(21b) 


We denote the vector v as the algebraic basis, in tha^ this variable represents the finite-dimensional 
approximation space with an equivalent number of < legrees-of-freedom; the proper functional basis 
corresponds to the images of = (1,0,0...),<0. 1.^ ),-,(0,...,0, 1) in v* through the transfor- 

mation Q, acting on the local baaes ais described by (14). 


3.2 Discrete Equations 

Armed with the variational forms of Section 2 Ji.nd the bases of Section 3.1, it is now a simple 
matter to construct the discrete equations. In parti cuIm, we note that our basis construction (21) 
allows us to express admissible elemental degrees-of- freedom v* in terms of v. This, in turn, permits 
us to construct the global discrete equations directly from local structure-preserving elemental 
equations, which is at the heart of the discretization-driven domain decomposition approach. 

We first construct the decoupled elemental matrices and inhomogeneity. 


blk{A^) = 


(Vhph„Vh<h,)‘=i 0 

0 (Vhph,,Vh.h,)*=* 


0 

0 


(Vhph„Vh.hy) 


k=K 


( 22 ) 
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mh 






V (M,,/)*=^ ) 

The A:th block of 6/A:(A*) represents the Neumann Laplace operator on the elemental domain fl*. 
We now recognize that not all elemental hihj are possible, and that not all hphg are admissible; 
indeed, the admissible degrees-of-freedom follow from the Q transformation of (21). We thus arrive, 
rather simply, at the fully discrete equations: 


Q^blk{A’‘)Qu = Q^blk{p). (23) 

We note that independent of the size of the mortar nullspace (Q right nullspace), (23) is solvable, A 
suflBcient condition for a unique mortar function is that Q'^Q be invertible; in the conforming cases 
Q^Q is simply the multiplicity of a node (that is, the number of elements in which it appears) . 

Equation (23) illustrates that the global Laplace operator can be thought of £ls a local operator 
“mortared” together by the Q'^jQ operations; indeed, the operator is the algebraic form of 
the standard direct stiffness procedure (here extended to nonconforming elements). In the imple- 
mentation of iterative procedures the are, of course, never explicitly formed, but rather are 

evaluated; diagrammatic evaluation of (direct stiffness summation) is shown in Figure 6 in terms 

of the diagram conventions defined in Table 1. The domain decomposition decoupling afforded by 
the implicit construction of the image basis through u* allows for efficient parallel implementation 
following the methods described in [11] for conforming techniques. 

Although the emphasis in the current paper is on the mortar discretization, the bases and 
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Figure 6: Direct Stiffness Summation of Residuals on Nonconforming Mesh 
evaluation procedure have been tailored to admit efficient iterative solution, and it is therefore 
appropriate to briefly indicate how the method is u led in conjunction with (for example) conjugate 
gradient iteration. To solve (23) we write 

r^ = Qnik{h-Q’Uk(A^)Qm So = 10 P'*) 

o™ = (r„,r„)/(s„.a’'W*:(A‘)03„) 


Um+1 — Mm 




im+1 = ^+1 

where m refers to iteration number, is the rrsidual, the search direction and (.,.) is the 
usual discrete inner product. All evaluations ar ? performed through the diagrams of Figures 4 
and 6. The blk{A^) operations are entirely local it the elemental level, with all transmission and 
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coupling through Q. The local A* cdculations are the standard conforming spectral element tensor 
product evaluations, as the mortar decoupling allows all local structure to remain intact despite 
global irregularity (e.g. non-propagating mesh refinement). 

4 The Stokes and Navier-Stokes Problems 

In this section we consider the extension of the nonconforming mortar method to the solution 
of the two-dimensional steady Stokes problem in a rectangularly-decomposable domain fl, 

-t/V^u-Vp = f (25) 

divu = 0, 

with homogeneous Dirichlet velocity boundary conditions u = 0 on dil. Here u is the velocity, p 
is the pressure, f is the forcing vector, and v is the kinematic viscosity. The associated vziriational 
problem is; Find (u,p) e (^fi(n)*, L*(n)) such that 

i/(Vu, Vw) - (p, divw) = (f,w), VwGH^(n)* (26) 

(g,divu) = 0, Vg€i*(n), 

where Z,*(fl) is the space of L* functions of zero mean. 

The discrete formulation of the problem consists of choosing two discrete approximation spaces, 
one for the velocity field and one for the pressure. It is shown in [3,22,24] for the conforming 
spectral element approximation that choosing both of these spaces to be polynomials of degree 
less than or equal to the same degree N leads to an ill-posed problem, in which spurious pressure 
modes arise [5,12]. The existence of such modes is in contradiction with the verification of the “inf- 
sup” condition [5]. As regards our nonconforming methods for the Stokes problem, our starting 
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point is the conforniing staggered mesh method defined in [23] for which the inf-sup condition is 
satisfied. The correct nonconforming extension is to ise the velocity space defined m (9a, 9b), 

and for the pressure the space Mh = {(f> e € PAf- 2 (n*)} associated with the conforming 

approximation. In essence, the fact the pressure is implies that it need not be modified when 

the constraints on the velocity are relaxed. 

With these spaces we arrive at the following noi conforming discretization: 

Find {uh,Ph) e i{Xhy,Mk) such that 

i/(Vuh, Vwh) - (pfc, divwh) = < f.wii >, Vwh G (27) 

(gfc,divuh) = 0, Vgfc G Mfc 

from which uniqueness, stability, and spectral erroi properties foUow firom the results of previous 
sections and [4], suitably modified within the Stokes context as described in [23,24]. (We note 
that, as elsewhere in this paper, we do not dwell on quadrature issues which are, by now, standard 
practice.) We then choose a basis for Mh 

i=i y=i 

where the are the N - 2th order Gauss-Legendre interpolants, that is, those polynomials of 

Fn-2 such that where an the W - 2 zeroes of ijv-2 [30]. We thus arrive 

at the discrete saddle problem 

Cf'blk{A^)Qn- Q^blkip^^p = Q^Wfc(f‘) (29) 

61jt(D*)Qu = 0. 

Here p is the algebraic basis for p analogous to u of (21b), Wfc(A*),Q^,g, WA:(f*) are defined as in 
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(18-22), and D is the gradient operator given by 

■ ( 9 ^-^ 9 ^ 0 . 0 

6/ifc(D*) = 

V,-.y € {0,...,JV}2, Vp,9 € {l,...,iV - 1}^ 

Extension to Navier-Stokes is straightforward given the lower-order nature of the convective terms. 

As in the pure elliptic discretization, (29) is amenable to iterative solution. We currently use 
a semi-implicit procedure for Navier-Stokes, in which the nonlinear terms are treated explicitly, 
and the Stokes subproblem is handled with a Uzawa nested iteration [20]; conforming multigrid 
techniques [28] are currently being extended to the nonconforming case. In addition to the staggered 
mesh Stokes treatment, elliptic-splitting methods appropriate for higher Reynolds number flows are 
also used [17,19]; these discretizations represent sequences of elliptic operations (23), and thus their 
extension to the nonconforming case follows from Sections 2 and 3. 

5 Numerical Examples 

In this section we illustrate various eispects of our nonconforming method by a number of exam- 
ples. The central point is the flexibility and ease-of-implementation afforded by a nonconforming 
approach based on a consistent, non-context-dependent matching. As our first example we consider 
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the Helmholtz problem 

-V*u + A*u = / cn n=]0,l[x]0.l[ (31) 

u = ca an 

where f is chosen such that the exact solution in n is given by u • In Figure 

7a we show a high-resolution conforming mesh h = {K = 16, N = •)] in Figure 7b we show a 
nonconforming mesh h = {K = 10,N = •),in which the local structure-preserving mesh refinement 
is illustrated. In Figure 8 we plot the error in the Xh norm of (10) for both solutions as a function 
of N {K fixed) for A = 50. This example demonstrates the rapid (here exponential) convergence 
of the spectral element approach, and the superior resolution properties of the nonconforming 
discretization, which achieves the same accuracy as he conforming approximation with significantly 
fewer degrees-of-freedom. 




Figure 7: High Resolution a) Conforming and b) Nonconforming Meshes for Helmholtz Problem 

As our second example we demonstrate the u ility of nonconforming methods in constructing 
appropriate meshes; we consider the labyrinth cha inel of Figure 9, in which the two meshes for the 
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O conforming (Fig. 7a) 

A nonconforming (Fig. 7b) 



Figure 9: Nonconforming Labyrinth Channel Mesh, h = {K = 22,N = 9) 

two sides of the channel are constructed “separately” , and subsequently merged by mortar. The 
boundary conditions are given as: a parabolic velocity profile at inflow, no slip on the channel walls, 
and outflow (constant pressure) at the exit. In Figure 10 we show streamlines for the steady Stokes 
flow calculated by the discretization (29) and the nested conjugate gradient Uzawa method; notable 
are the continuity at element boundaries and the lack of spurious pressure modes. The mesh in 
Figure 9 can be thought of as one instance of a sliding channel calculation; nonconforming methods, 
with appropriate extension (as in Figure 1), should prove to be powerful techniques for moving 
boundary problems when used in conjunction with arbitrary-Lagrangian-Eulerian techniques [15]. 



Figure 10: Stokes Solution for the Labyrinth Channel 

Lastly, we consider a moderate Reynolds number flow past a wedge [17,26] in a channel; the 
utility of the nonconforming methods in generatin g an appropriately refined mesh is apparent in 
the mesh shown in Figure 11. Note that we relax here the constraint, introduced for simplicity 
of exposition in previous sections, that the elemen jS be rectangular; treatment of general curved 
elements represents a simple extension of the methods described in Sections 2-4. In Figure 12 we 
show the short time solution of the startup vortex n-3ar the tip of the wedge, for a Reynolds number 
R = 500 at a time r = ^ = .085 on the mesh h = = (if = 16,iV = 9) of Figure 11. We prescribe 
a slug velocity profile at inflow, no slip boundary conditions on the walls, and outflow (constant 
pressure) at the exit. Here R = where V is the channel average velocity, H the channel width, 
and 1 / is the kinematic viscosity, and t is time. The high resolution in the vicinity of the wedge 
allows for a detailed description of the startup vor ex. 
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